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We present a global fit to the structure function F2 measured in lepton-proton experi- 
ments at small values of Bjorken-a;, x < 0.01, for all experimentally available values of 
Q 2 , 0.045 GeV 2 < Q 2 < 800 GeV 2 , using the Balitsky-Kovchegov equation including 
running coupling corrections. Using our fits to F2, we reproduce available data for Fl 
and perform predictions, parameter-free and completely driven by small-x evolution, 
to the kinematic range relevant for the LHeC. 

Deep Inelastic lepton-hadron Scattering (DIS) experiments constitute one of the main 
sources of information about the partonic structure of hadrons, as well as a testing ground 
of pQCD techniques. The DIS program carried out in HERA has explored an unprece- 
dented kinematic regime in the DIS variables Q 2 and Bjorken-x. At small values of x, 
the emission of additional soft gluons in the hadron wavcfunction is enhanced by factors of 
a s In l/x ~ 0(1), resulting in a fast growth of the gluon density of the hadron with decreas- 
ing Bjorken-x. Such growth has been observed in data and, at small enough values of x, or, 
equivalently, for large gluon densities, it is expected to be tamed by repulsive gluon-gluon 
self-interactions. The latter mechanism is known as saturation of the gluon distribution. It 
is non-linear phenomenon, closely related to unitarity of the theory, which is not included 
in the more standard linear approaches to DIS based on the DGLAP or BFKL pQCD evo- 
lution equations. Here we present an analysis of available DIS data at small-x based on 
the non-linear Balitsky-Kovchegov (BK) equation including the recently calculated running 
coupling corrections to the evolution kernel [TJ [3J [3J 0] . 

The starting point of our analysis is the dipole model of DIS, which stems from the 
kt -factorization theorems in the high-energy limit. It states that, at small-x, the virtual 
photon-proton cross section can be written as a convolution of the light-cone wavefunction 
for a virtual photon to split into a quark-antiquark dipole, I^t,!,! 2 , times the cross section 
for dipole-proton interaction. For a hadron homogeneous in the transverse plane, one has: 

a^l(x,Q 2 )^a Q f dzdr\* T , L (z,Q 2 ,r)\ 2 Af(r,x) , (1) 
Jo 

where the subscripts T, L stand for the polarization of the virtual photon (transverse or lon- 
gitudinal), r is the dipole transverse size and Q 2 the photon virtuality. Explicit expressions 
for ^t,l hi the case of any number of active flavors can be found in [5] . In our analysis we 
consider only the three flavor case, i.e, with no charm quark. All the information about the 
strong interactions, and also all the x dependence of the absorption cross section, is encoded 
in the imaginary part of the dipole scattering amplitude, Af(r, x) (related to the dipole cross 
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section through the optical theorem). The normalization a® in Eq. (JTJ) originates from inte- 
gration over the transverse plane, trivial in the case of a cylindrical proton considered here, 
and will be one of the free parameters in the fit. The inclusive DIS structure function is 
then given by 

F 2 (x,Q 2 ) = -4 (°t + <jl) , (2) 

4 ir^a, 



em 
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whereas the longitudinal structure function Fl(x,Q ) is obtained by just considering the 
longitudinal component. 

General arguments in QCD indicate that the dipole scattering amplitude is small for 
small dipole sizes r (color transparency), whereas large dipoles suffer a large absorption. 
Phenomenological works in the context of saturation parametrize the transition between 
these two regimes, dilute and saturated, in terms of a i-dependent saturation scale Q 2 s (x), 
such that N(r,x) ~ r 2 for r <C 1/Q s (x), whereas Af(r,x) ~ 1 for r > 1/Q S - The details 
of the paramctrization depend on the physical input considered in each case. The leading 
cc-dependence of the structure functions is thus encoded in that of the saturation scale, com- 
monly assumed to have a power-law behavior: Q 2 (x) ~ (xo/x) x GeV 2 , with the evolution 
speed A being a free parameter fitted to experimental data. The results from different groups, 
albeit using slightly different parametrizations of the dipole amplitude, report A ~ 0.2-^0.3. 

In this work, rather than resorting to phenomenological models, we calculate the x- 
dependence of the dipole amplitude in terms of the pQCD non-linear BK evolution equation 
including the recently calculated running coupling corrections. Importantly, the original, 
leading- log (LL) in a s lnl/x BK equation yields an evolution speed A ~ 4.8 a s , with a s the 
bare coupling \§\ . Such value is too large to attempt a good description of data. Importantly, 
the running coupling corrections to the LL kernel obtained through resummation of a s Nf 
terms bring the evolution speed down to values compatible with those extracted from data, 
among other interesting dynamical effects [4] . The BK equation reads 

cW(r, Y) 
d \tixq/x 

x [JV (n , x) + N(r 2 , x) - N{r, x) - M{n , x) N(r 2 ,x)] , (3) 

where = r — ri and xq is the starting point for the evolution. In our case xq = 10~ 2 
is the largest experimental value of x included in the fit. The precise shape of the running 
coupling kernel K depends on the scheme chosen to single out the ultra-violet divergent 
contributions from the finite ones that originate after the resummation of quark loops. It 
was found in [4] that the prescription suggested by Balitsky in [2j minimizes the role of the 
sublcading contributions, and is therefore better suited for phenomenological applications. 
The corresponding kernel is 



J drx K(r,r 1 ,r 2 ) 



N c a s (r 2 ) 
K(r,n,r 2 )= ^ 



1 (ccsirl) A 1 (a s {rl) 1 



1 ot s {rl) J r 2 \a s (r 2 ) 



(4) 



Details about the numerical method used to solve the BK equation can be found in [7]. 
We evaluate the coupling according to the standard one-loop QCD expression. In order to 
regularize it in the infrared, we freeze it to a constant value for r > rf r , with ot s (r 2 r ) = 0.7. 
Thus 

a s (r 2 ) = for r <r fr , (5) 
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Initial condition 


cr (mb) 


Qlo (GeV 2 ) 


C 2 


7 


X 2 /d.o.f. 


GBW 


31.59 


0.24 


5.3 


1 (fixed) 


916.3/844=1.086 


MV 


32.77 


0.15 


6.5 


1.13 


906.0/843=1.075 



Table 1: Values of the fitting parameters from the fit to available f^x, Q 2 ) data with 
x < 10~ 2 and for all available values of Q 2 , 0.045 GeV 2 < Q 2 < 800 GeV 2 . 




Figure 1 : Comparison between a selection of experimental data and the results from the fit 
for F2(x, Q 2 )- Solid red lines correspond to GBW i.e., and dotted blue ones to MV i.e. 



and a s (r 2 ) = 0.7 for r > rf r . The parameter C in Eq. © is a free parameter in the fit. We 
take Aqcd = 0.241 GeV. A final ingredient in our analyses is the initial condition for the 
evolution, i.e. the specific form of the dipole amplitude at x = 10~ 2 . We considered two 
different initial conditions, GBW: 



A" 



GBW 



(r, x = 10 2 ) = 1 - exp 



Q 



s0 



and MV: 



N MV {r,x = 10~ 2 ) = 1 - exp 



?s 



hi 



1 



+ e 



(6) 



(7) 



, r Aqcd 

where Q 2 is the initial saturation scale and 7 is another free parameter which controls the 
small r behavior of the dipole amplitude. For the fit with GBW i.e. 7 is kept fixed to 1. 

All in all, there are four (three when 7 is fixed) free parameters in the fit: The ini- 
tial saturation scale Q s q, the overall normalization ctq, the infrared parameter C and the 
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anomalous dimension 7. With this setup we obtain an excellent fit to available experimen- 
tal data for the inclusive DIS structure function F2(x,Q 2 ) with x < 10~ 2 and all existing 
values for Q 2 , 0.045 < Q 2 < 800 GeV 2 , which adds up to 847 data points. In order to 
smoothly go to photoproduction, we have used the following redefinition of the Bjorken 
variable: x = x (1 + 4m 2 /Q 2 ), with to/ = 0.14 GeV for the three light flavors considered in 
this work. The experimental data included in the fit has been obtained by the HI and ZEUS 
Collaborations (HERA) and by the E665 (Fermilab) and NMC (CERN) Collaboration^]. 
The parameters from the fit are presented in Table 1. We have checked that the quality of 
the fit and the values of the parameters are stable under the restriction of the data range 
to the region Q 2 < 50 GeV 2 (which leaves 703 data points for the fit). A comparison be- 
tween a selection of experimental data and the results from the fit for F2(x,Q 2 ) is shown 
in Fig. [1] The two initial conditions GBW and MV offer equally good fits to data. With 
the parameters values yielded by the fit, we obtain a good description of the longitudinal 
structure function, as shown in Fig. [2j 



data: H1 (PLB665, 139; x-averaged) 
solid: GBW initial conditions 
dotted: MV initial conditions 
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data: ZEUS (0904.1092; unconstr. errors) 




solid: GBW initial conditions 
dotted: MV initial conditions 
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Figure 2: Comparison between experimental data from the HI (left) and ZEUS (right) and 
the predictions of our model for Fl(x,Q 2 ). The theoretical results have been computed at 
the same {x) as the experimental ones. 

The analyses of available data, including this work, offer compelling evidence for the 
presence of gluon saturation effects. However, a clear distinction from alternative (linear) 
approaches cannot yet be established unambiguously. Indeed, a good description of all 
data for Q 2 > 1 can be achieved using DGLAP fits, even though the very application 
of DGLAP at such low Q 2 is questionable. On the other hand, the obscrvables like Fl 
which are more sensitive to the gluon distribution are poorly determined experimentally. 
The proposed LHeC [5] may have the potential to disentangle different physical scenarios 
through a deeper exploration of the small- a; kinematic region. The question rises of which 
observable is best suited to pin down the saturation of the proton gluon distribution. A 
partial answer is given in Fig. O courtesy of Juan Rojo [TT]. It presents a NNPDF1.2 
DGLAP fit [jJJ [TU] to pseudodata for proton F 2 and F ^ in the kinematic regime relevant for 
the LHeC (down to x = 10~ 6 approx.). The pseudodata are generated by extrapolation of 
the fit to F2 presented here using BK evolution, and thus include saturation effects. While 

a A detailed list of references to the publications where data is presented can be found in [7] 
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Figure 3: NNPDF 1.2 DGLAP fits (green stars) to pseudodata (red bars) for proton F2 
(left plot) and F^ (right plot) generated by extrapolation of our fits down to x = 10~ 6 for 
Q 2 = 2,5,10 and 20 GeV 2 . 

F2 can be well fitted by DGLAP, the best DGLAP fit to pseudotata underestimate Fl at 
small Q 2 and overshoots it at the largest Q 2 considered, yielding a large x 2 /d.o.f.. One 
concludes that a precise experimental determination of Fl at the LHeC over a large enough 
Q range might suffice to pin down the kinematic region where departure between DGLAP 
and non-linear (or, maybe, linear resummcd small- a; evolution) takes place, while F2 does 
not offer such discrimination power. 
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